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ABSTRACT 

The maser disk around the massive black hole (MBH) in active galaxy NGC 4258 
exhibits an 0(10°) warp on the 0(0.1 pc) scale. The physics driving the warp are still 
debated. Suggested mechanisms include torquing by relativistic frame dragging or by 
radiation pressure. We propose here a new warping mechanism: resonant torquing of 
the disk by stars in the dense cusp around the MBH. We show that resonant torquing 
can induce such a warp over a wide range of observed and deduced physical parame- 
ters of the maser disk. 

Subject headings: galaxies: nuclei — stellar dynamics — galaxies: individual (NGC 
4258) — accretion disks — masers 



Introduction 



There is compelling evidence that massive black holes (MBHs) of rnass 10^ Mc^ < M. < 
few X 10^ Mc^ exist in the centers of most galaxies ( Ferrarese & Merritt 2000 : Gebhardt et alj|2003 : 
Shields et alb oosk T he MB Hs acquired most of their mass by efficient luminous accretion (|Soltanl 
19821 : Iyu & T remai nell2002h. which implie s accretion from an optically thick, geometrically thin 



accretion disk (jshakura & Sunyaev 1973). On larger spatial sca l es, such disks may fragment 



and form stars o n disk-like orbits (IKolykhalov & Svunvaevlll980l: IShlosman & Begelman 



1989; 



Goodmanll2003h . as o bserved in the Gala ctic Center (iBeloborodov & LevinI 120041 : iPaumard et al. 
20061) and elsewhere (|Pizzella et al.ll2002l) Gas disks, often m arginally stable and with properties 
similar to star disk progenitors (IMilosavljevic & Loebll2004h . are dete cted in the radi o by H2O 
maser emissio n (for a list of kno\yn cir cum-nuclear maser disks, see iMaloneyl l2002l : for more 
candidates, see lBraatz & Gugliuccill2008|) . 



The interest in processes that affect accretion disks stems from the possible implications for 
MBH growth, for the luminosity evolution of Active Galactic Nuclei (AGN) and for nuclear star 
formation. Conversely, the disk properties may provide information about dynamical processes 
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that operate very close to the MBH, and may be relevant for other phe nomena there, such as tidal 
interactions with the MBH or the emission of gravitational waves (e.g. lAlexandeiil2005L |2007|) . 



Maser-emitting nuclear accretion disks are unique, clean probes of the near environment of 
MBHs. The first discovered and be st studied mase r disk, NGC4258, is also the thinnest and 
most Keplerian (to better than 1%, iMaloneyl l2002h nuclear disk yet discovered. Radio H2O 
maser observations of the disk morphology, velocities and acc elerations allow accura te measure- 
ments of the MBH mass and the distance to the host galaxy (IHermstein et al.lll996h . Together 
with opti cal and X-ray observation s, these provide estimates the mass accretion rate through 
the disk (INeufeld & Malonevlll995l). T he NGC 4258 maser disk shows a clear 0(10°) warp on 
the fewxO.lpc scale (|Herrnstein et al.l Il996|) . whose orig in is still a matter of debate. Warps 
are possibly also obser ved in the maser disks of Circinus (jGreenhill et al.ll2003|) and NGC 3393 
JKondratko et allboosk 



Several mechanisms were proposed to explain the warped disk in NGC 4258. The lBardeen & Petterson 



(| 19751) (BP) effect (torquing of a viscous disk by the General Relativistic frame-dragging) can warp 
a disk that is initially misaligned with the MBH spin by dragging its inne r regions to the MBH equa- 
torial plane ([Lodato & Pringlelbood: ICaproni et aklboovl : iMartinlboOSI) . Instability of a thin disk 
to rad iation pressure from the central sou rce can result in the amplification of a small pre-existing 



warp (|Pringlell 19961 : Malonev et al.ll 19961). The disk's self-gravity can support equilibrium warped 
disk configurations (|Ulubay-Siddiki et al.ll2007|) . if these are excited by other means. One feature 
shared by these mechanisms is that they require some initial asymmetry in the disk or disk/MBH 
alignment for the warp to develop. 

Here we consider the implications of the dense stellar cusp around the gas disk. Various dy- 
namical scenarios predict the formation of steep stellar density cusps n(r) oc (r/r/,,)""^ (0.5 < 
7 < 2.5) within the MBH radius of influence (e.g. iBahcall & Wolf] Il977|) . The empirical 
M,/a correlation implies that the cusp density scales as Mr^^^, so that the relativel y low-mass 
MBH s, where maser disks are fo und, are surrounded by a v ery dense stellar cusp ([Alexander 
20071). Strong mass segregation ([Alexander & HopmanI 120091) further concentrates stellar-mass 
black holes (SBHs) to the center. The stars in the cusp move in a nearly spherical potential dom- 
inated by the MBH, and therefore orbit on nearly fixed planes. This syr nmetry leads to the rapid 
torquing of any test particle by the p rocess of resonant relaxation (RR) ([Rauch & Tremainelll996 
(RT96); iHopman & Alexandeiil2006h . In particular, gas streams on circular orbits undergo "vec- 
tor RR", which changes the orbital plane, but not their eccentricity. RR thus induces warps by 
exchanging angular momentum between the stars and the disk. 

This letter is organized as follows. RR dynamics are summarized in ^ An analytical model 
for RR torquing of a disk is described in ^ and applied to the NGC 4258 maser in ^|4l Our results 
are presented in ^|5]and summarized in ^ 
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2. Resonant relaxation dynamics 

RR is a rapid relaxation mechanism of the angular momentum L, which operates in near- 
symmetric potentials that suppress the evolution of the stellar orbits (e.g. fixed ellipses in a Kepler 
potential, or planar rosettes in a spherical potential). In such systems, the residual torque on a test 
mass by the orbit- averaged mass distributions of A^^ stars of mass Af^ within distance r from the 
center, |T| ~ N^^GM^/r, remains constant on timescales shorter than the coherence time to, as 
long as perturbations due to deviations from the perfect symmetry remain small. The change in L 



then grows coherently as oc t on timescales P <^t <^to, where P = I-k^t^ jG^M, + N-^M^) is 
the circular orbital period. When t ^ t^, the large accumulated change over the coherence time, 
(AL)o = I T|to, becomes the step size of a rapid random walk, (AL/Lc)(t) = [(AL)o/Lc] yJtJtQ = 
\/t/T^, where Lc = \JG{M, + N^M^)r is the circular specific angular momentum for the test 
mass' orbital energy, and Trr is thereby defined as the RR timescale, when \AL\/Lc = 1. The 
slower the loss of coherence ("quenching"), the more efficient is RR (shorter Trr). RR can be 
orders of magnitude faster than non-coherent 2-body relaxation. 

In a near-Keplerian potential the orbits are nearly fixed ellipses, and the coherence time is set 
by the faster of General Relativistic (GR) precession or precession by the potential of the enclosed 
stellar mass. In this case RR can change both the direction and magnitu de of L ("scalar RR"). Ho\y - 



ever, the torques along L fall to zero for a test mass on a circular orbit (IGiirkan & Hopmanll2007l) . 
and so circular orbits, such as those of gas streams in an accretion disk, remain circular (ATy = 0), 
but their orientation evolves rapidly (AT^ ^ 0, "vector RR") until RR itself randomizes the orbital 
planes to a degree where coherence is lost ("self-quenching"). Vector RR is not quenched by pre- 
cession (AT^ 7^ between orbit-averaged rosettes), and so it can operate both very near the MBH, 
where GR precession is fast, and far from the MBH where the potential is spherical but no longer 
Keplerian. The coherence time for the self-quenched vector RR is to = ^o-^c/|T| = Aq-/NIP/ /i, 
where ^ = N^^M^/ (M, + N^Mi^) and where Aq is an order unity factor. On timescales t<^tQ, 

' - P.^V^A (iV.M.«M.). (1) 



where is an order-unity factor (recent A^-body simulations indicate that Py ~ 1.8, Eilon et al 



20091). The vector RR timescale is then 



A^y^ „ 1 fM,\ p 



(3^Ao /i (3^Ao \M^J y/N^ 

For the purpose of numeric evaluation, we assume below P^Aq = 2. Note that for vector RR, 

^vRR~^0- 
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A change ofw = | AL|/Lc in the angular momentum of a gas ring in an accretion disk ("warp 
factor" w), expressed by the change to the ring normal An= {9, ip), corresponds to a warp inclina- 
tion angle 6 via = 2(1 — cos 6), For w <^l,w ^ 6 rad. By definition, 6 = 60° {w = 1) over the 
vector RR timescale Tvrr. RR warping by a factor w occurs on a timescale shorter than the 
coherence time {w oc t), and so the time required to develop such a small local warp inclination is 

t(w)~Wto~W^vRR- 

The residual RR torque T(r) changes gradually with r, as the uncorrelated stars that ef- 
fectively contribute to it vary. We estimate the RR spatial coherence factor = r2/ri as the 
distance ratio over which the enclosed stellar numbers satisfy 2y^N{ri) = N{r2), which trans- 
lates to fc = 22/(3-7). The observed NGC4258 maser disk extends over R2/R1 = 2, and the 
assumed 7 = 1.5 stellar cusp there corresponds to /c ~ 2.5 which implies that the warping 
torques, and thus Ani and An2, are approximately uncorrelated (section |4]). The total observed 
warp angle cu = cos^^(Ani ■ An2) is due to the difference between the local warp angles at the 
disk's inner and outer edges. The relative magnitude of the warp factors is related by Eq. ([T]), 
U = W2/W1 = {R2/Ri)-^'^ (U = 0.6 for 7 = 1.5). In the limit of small warp inclination an- 
gles, as is the case here, f^j ~ 6'2/6'i and uj ~ 6'ia/I — 2/^ cos(v9i — 922) + fw- The maser disk in 
NGC 4258 dis plays a = 8° relative warp between the limits of the masing region at Ri = 0.14 pc. 



i?2 = 0.28 pc ( Herrnstein et al. 19961) . which after averaging over uncorrelated ipi and (p2, implies 



that locally, Oi ~ 10° and 62 ~ 6°. 



3. Accretion disk warping by resonant relaxation 



The gas in a thin accretion disk flows slowly into the MBH on nearly circular orbits. Over 
time, the orbit associated with a ring element of the disk shrinks, and the torques on it change. RR 
can significantly affect the ring only if the typical RR timescale, T^rr is shorter than the radial 
inflow timescale U ~ R^ /y, which is set by the disk's tangential in-plane viscosity coefficient, 
z/, where R is the radius along the disk's mid-plane. The presence of a warp also induces a ver- 
tical shear vis cosity vo, which resists the external to rques by diffusing the warp on a timescale 
tyj = R^ /i'2 (|Pringlelll992l: iLodato & Pringlell2007h . Effective warping thus also requires that 
^vRR be shorter than the warp diffusion timescale. In addition to these timescale constraints, the 
torquing stars must carry enough residual angular momentum to exchange with the disk and warp 
it. Therefore, both a timescale condition ( §3.11) and an angular momentum condition ( §3.21) must 
be satisfied for RR torquing to occur. 
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3.1. RR timescale condition 



The vector RR timescale (Eq. ^ estimated above for stars enclosed within radius r and acting 
on a mass orbiting at the same typical radius r, can be easily generalized to a mass (disk ring 
element) at a different radius R, by noting that T^rr scales with the torque as oc |T|^^, and that 
the typical torque on the ring element is decreased by a factor r/R when R>r and by a factor R/r 
when R<r, 



TvRR(r, R) = T.RR(r) (r//?)^® oc r^^+^/s ^ (3) 

where 6 = sgn(r — R), and the last proportionality holds in the Keplerian limit. In order to 
efficiently warp the disk by a factor w, T^rr must be shorter by at least some factor < 1 than the 
faster of the two viscous processes, Uw = min(tr, tw). 



ettrM>wT^RR{r,R). (4) 

An equality in Eq. (01) defines an upper limit r'f''' for a given R<r and a lower limit r[^^ for R>r 
(assuming < 7 < 4; recall that 0.5 < 7 < 2.5 for realistic cusp models) 



,(e) 



Iso 



-M, {et/w)trUR) (R^ 

M. Pi,o Ulsc 



20 



1/(26+7/2) 



(5) 



where A^iso is the formally extrapolated number of stars within the last stable circular orbit riso, and 
Piso is the formal Keplerian orbital period there. When there exists a solution with r\~^ < r'f'\ then 
stars in the volume rj^-* < r < r^^-* induce fast enough vector RR on the disk at R to warp it by a 
factor w. 



Note that the timescale condition implies that the RR torque dominates over the viscous 
torques, since |Trr| ~ Lc/TvRR>Lc/t™~max(|Tr|, |T^|). 



3.2. RR angular momentum condition 

Efficient torquing of the disk mass in a ring (i?, R+/S.R) by stars in a spherical shell (r, r + Ar) 
requires that the disk's angular momentum be smaller by at least some factor < 1 than the 
random residual angular momentum carried by the stars, assumed to be isotropically distributed 
(this a conservative criterion, since it neglects the possibility of warping by transfer of angular 
momentum from the disk to the stars). The total angular momentum in both the disk and the 
stellar cusp is dominated by the large scales (e.g. a disk with surface density (x R^ ( ^4.21 ) and 
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temperature that falls with i? has s > -3/2 (Eq. ^ and LaAMa oc i^^/^+^Ai?, where is the 
disk mass enclosed within R; a stellar cusp with 7 < 2 ( ^4.11) has L^A(-\/iV^M^) cxr i-T/2Ar). We 
therefore approximate the angular momentum criterion by requiring that the angular momentum 
in the disk within R be smaller than the residual angular momentum in the stars within r, 



wMa{R)L,{R) < eL^/fUr)M,L,{r) cc {r/n,,f-^'^ . 
An equality in Eq. ^ defines a lower limit ri (assuming 7 < 4) 

-n 1/(2-7/2) 



r\so 



MdiR) 



R_ 

r\so 



(6) 



(7) 



When there exists a solution such that ^ = maxfr.^ \ 



ri) < r 



(+) 



then the stars in the volume 



.(-) 



< r < r 



(+) 



carry enough residual angular momentum to significantly torque the disk enclosed 
within R, on a short enough timescale. 

Note that the angular momentum condition implies that the RR torque dominates over the 
disk's self-gravity, since the self-gravity torques \Td\ ^GMd/r<G^/NlMi,/rr^ |Trr|- 



4. NGC 4258 model 

In order to apply the RR torquing mechanism to the maser disk in NGC 4258, it is necessary 
to specify the stellar cusp density profile, the disk's viscous timescales, trw{R) (Eq. S]), and its 
enclosed mass Md{R) (Eq. |6l). We adopt here simple, observationally motivated models of the 
cusp and disk. 



4.1. Galactic nucleus model 



NGC 425 8 is a spiral galaxy at a distance of 7.2 ± 0.3 Mpc wi th a M. ~ 3.7 x 10'^ 



central MBH ( Herrnstei n et al.lll999l). T he empirical M./cr relation dFerrarese & MerrittI I2OO0I : 
Gebhardt et allEoOSr : iShields et al.l boosh implies that a ~ 150 km s"^ and the MBH's radius of 
dynamical influence is rh = GM,/a'^ ~ 7 pc, where the enclosed stellar mass is UhM,, with Uh 
0(1) . On the distances spanned by the maser spots, Ri = 0.14 to R2 = 0.28 pc (IHerrnstein et al. 
1996h . main-sequence stars dominate the population with a relatively flat 7 ~ 1.5 power-law density 



profile, typical of the low-ma ss component in a mass- segregated population ( Bahcall & WoHlll977 : 
Alexander & Hopmanll2009l) . We model the central cusp of NGC 4258 as a power-law cusp, A^*(< 
r) = /i/j(M,/Mj,)(r/r/i)^~'^, with 7 = 1.5, r/^ = 7 pc, yUh = 2 (the formal value for an MBH-less 
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singular isothermal distribution) and = 1 Mq stars. The corresponding vector RR timescale is 
TvRR~fewx 10'' yr across the disk. 



4.2. Maser disk model 



Following ICaproni et al.l (|2007|) . we assume a power-law surface density profile for the disk, 
S(i?) = T.i{R/Riy. We assume the structure equations of a stationary, georn etrically thin, opti- 
cally thick Keplerian disk around a MBH of mass M, (e.g. Frank et al.ll2002|) . These can be ex- 
pressed in terms of the mass accretion rate M, the dimensionless viscosity parameter a = u /cgH, 
the surface density power-law index s and p{Ri), the mid-plane mass density of H2 at Ri (the disk 
mass is assumed to be mostly in H2). The normalization Si at Ri is fixed by 



S^ = Mp2 /2an 



(8) 



where f2 = ^JGM,/ R^ is the Keplerian frequency, H is the disk scale-height H, and is isother- 
mal sound speed c^. 



H^ = M/2Tram, c^ = Mn/27raS . 



(9) 



For an ideal gas = ^JuTJ^ with T the gas temperature and [i the mean molecular mass {\i = 2mp 
assumed). The radial inflow and the warp diffusion viscous timescales are then 



R^ / oiCgH , tyjr^ R^ / a2CsH , 

where ^2 = max[/(a), a2m»^] with f(a)c^2(l + 7a'^) /\a(A+a'^)] Jogilvielll999| ) and 02 max ' 
with some uncertainty (ILodato & Pringlell2007|) . 



(10) 

'3-4, 



Caproni et al.l (l2007h combine observational and theoretical constraints and deduce that the 
physical parameters of the NGC4258 maser disk lie in the range 10^® < eM < lO~'^M0yr^\ 
where e is the radiative ac cretion efficiency (but see higher estimate M > 10^'^ Mq yr~^ for ADAF 



model, iLasota et al.lll996|) . and 0.03 < a < 0.2. They consider a range of possible disk profiles, 
— 2 < s < 0. Efficient production of H2O maser emission constrains the disk's physical parameters: 
the H2 density, 10^ cm"^ < < 10^^ cm"^, t he gas temperatu re, 400 K < T < 1000 K, and the 
gas pressure 10^*^ K cm"^ <p/k< 10^^ K cm" -3 ([Maloney"2002'). In addition, the disk aspect ratio 
is constrained by observations to H/R < 0.002 (Moran .20081) . The disk is gravitationally stable 
only if Md/M, < H/R. 
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5. Results 

Figure ([T]) shows the region in the disk's {M, a, s,nH2) parameter space where masing is 
possible with the aspect ratio and disk mass constraints ( §4.21) . and where both the timescale (Eq. 
IHwith et = 1, 02 max = 3), and angular momentum conditions (Eq. [6] with = 1) are satisfied for 
9 > 10° between 0.14 and 0.28 pc, so that RR torquing can warp the disk by at least 8°, on average. 
We find that such an RR-induced warp is possible over a wide range of the observed and deduced 
physical parameters of the maser disk. RR warping is less efficient for disks with a very steep 
density profile, high density or a high mass accretion rate. We conclude that the RR mechanism 
can drive a moderate warp in the NGC 4258 maser disk, but cannot substantially perturb or disrupt 
it. We also find that the highest acc retion rate con sistent with masing is M < 10^^ Mq yr~^, which 



agrees with X-ray based estimates (IWatsonll2002|) . 



6. Summary 



We have shown that the torques exerted by the Poisson fluctuations in the distribution of the 
stars around the MBH in NGC 4258 can transfer momentum from the stars to the maser disk and 
excite the warp observed in the disk on a timescale w Tvrr ~ 0(10^ yr), more efficiently than 
the BP effect (^rp >feAy x 10^ yr. ICaproni et al.ll2007|) or radiation pressure ( IGammie et al.ll 19991 : 
see ^ Caproni et al. 2007 and references therein). RR is inherent to the discreteness of the stellar 
system, and thus does not depend on particular initial conditions for the disk, such as the initial 
warp required for radiation pressure warping, or the initial misalignment between the disk and the 
MBH spin axis required by the BP effect. RR-induced warps are transient and vary on the RR 
timescale T^rr ~ O ( 10 yr) . 

The NGC 4258 maser disk stands out in its well-defined morphology, Keplerian dynamics and 
the available detailed high-quality observations. These make it ideal for testing the RR torquing 
model. However, vector RR operates also in non-Keplerian spherical potentials, and can warp any 
disk configuration, as long as its mass is small enough and the viscous timescales are long enough. 
It therefore likely that RR torquing is relevant in other maser disk systems as well. 

Generally, RR torquing of accretion disks is also ex pected on much smaller scale s, driven by 
the strongly segregated cusp of stellar mass black holes (I Alexander & HopmanI 120091) . RR could 
thus potentially affect the accretion rate and direction, and thereby MBH mass and spin evolution. 



We thank S. Tremaine, J. -P. Lasota and the anonymous referee for helpful discussions and 
comments. TA acknowledges support by ISF grant 928/06 and ERC Starting Grant 202996. 
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Fig. 1. — RR disk warping in NGC 4258 disk parameter space (a, M, s) . Left: Range of surface 
density slopes s with riHa (Ri) = 10^ cm^'^. Parameter space regions where RR warping by > 8° 
between 0.14pc and 0.28 pc can be achieved as observed, are marked by open circles, regions 
where RR is not efficient enough are marked by filled circles, and regions where the physical 
conditions in the disk are incompatible with H2O maser production ( ^4.21 ) are marked by dots. 
Right: The same, for s = -0.75 (SS73 a-disk model) and for rinARi) = 3 x 10^ 10^ and 3 x 
10^ cm~^. 



